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FOREWORD 


The  azimuth  and  elevation  of  the  sun  are  routinely  determined 
by  a  navigator  in  the  air  and  on  the  sea  to  find  his  vehicle's  di¬ 
rection  and  line  of  position  (LOP).  A  sun  LOP,  when  used  in  conjunc¬ 
tion  with  another  LOP  obtained  from  other  celestial  bodies  (moon  and 
Venus)  or  radio  navigational  aids,  enables  the  navigator  to  locate 
his  vehicle's  position  (latitude  and  longitude). 

The  method  described  in  this  report  to  find  the  sun's  azimuth 
and  elevation  angles  is  based  on  a  procedure  which  the  author  learned 
during  undergraduate  Navigator  Training  at  Mather  Air  Force  Base,  CA. 
It  is  such  an  easy  method  that  in  a  minute's  time  one  can  obtain  the 
azimuth  and  elevation  angle  of  the  sun  by  using  two  tables:  The  Air 
Nautical  Almanac  and  the  H.O.  249  Sight  Reduction  Table. 

Rather  than  manually  extracting  data  from  the  tables,  the 
necessary  celestial  mechanics  and  solar  motion  data  have  been  pro¬ 
grammed  into  a  computer  program,  SAZEL.  One  need  only  know  the  time 
of  day,  position,  and  date  to  obtain  positional  data  for  the  sun. 
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INTRODUCTION 


To  obtain  data  on  rotational  motions  of  projectiles  in  free  flight, 
a  device  called  a  "yaw  sonde"  is  used.  The  working  details  of  the  yaw 
sonde  are  described  in  references  1  through  4.  Attitude  data  are  ob¬ 
tained  during  the  flight  of  a  missile  by  making  repeated  measurements 
of  the  angle  between  the  longitudinal  axis  of  the  missile  and  the  sun. 

Any  change  in  this  angle  can  be  approximated  by  a  change  in  attitude 
of  the  missile,  since  the  sun  moves  infinitesimally  during  a  projectile's 
flight  time. 

In  order  for  the  yaw  sonde  to  obtain  measurements,  it  is  necessary 
that  the  sun  be  within  the  field  of  view  scanned  by  the  yaw  sonde  during 
the  flight  of  the  missile.  The  interval  between  the  first  and  last  times 
the  sun  appears  in  the  field  of  view  for  the  interesting  parts  of  the 
trajectory  on  a  given  day  and  at  a  given  location  is  the  "window".  To 
determine  the  window,  the  sun  position  (consisting  of  azimuth  and  ele¬ 
vation)  must  be  known  or  computed  from  other  sun  data.  References  5 
through  7  present  various  computer  methods  of  determining  the  window. 
Common  to  every  method,  however,  is  a  requirement  for  selected  input 
data  relative  to  the  sun,  which  is  obtained  from  references  8  and  9. 

This  report  describes  various  methods  of  determining  the  sun's 
azimuth  and  elevation  (for  any  year  from  1965  through  2000)  and,  by 
the  use  of  the  computer  program  SAZEL,  eliminates  the  requirement  for 
input  data  obtained  from  references  8  and  9. 

DISCUSSION 


Sun  Features 

The  sun  exhibits  certain  features  that  contribute  to  its  use  in 
the  yaw  sonde  method  of  measuring  missile  attitude: 

1.  Its  position  can  be  determined  accurately  at  any  instant 

in  time. 

2.  A  change  in  its  position,  in  relation  to  the  earth,  is 
of  so  short  a  duration  as  to  be  negligible  during  the  flight  of  most 
missiles. 

3.  It  provides  ample  energy  for  the  functioning  of  the  yaw 

sonde. 
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Sun  Acquisition  by  a  Yaw  Sonde 

The  signals  produced  by  the  yaw  sonde  can  be  reduced  to  the  spin 
rate  of  the  projectile  and  the  angle,  y,  between  the  sun  vector  and  the 
longitudinal  axis  of  the  projectile;  or  the  angle,  ct,  between  the  sun 
vector  and  a  norroal-to-the  projectile  axis  (Tig  2).  The  convention 
adopted  here  is  that  a  is  positive  if  the  sun  is  toward  the  nose  of 
the  projectile  and 

a  =  90°  -  y  Q.) 

The  ARRADCOM  yaw  sonde  is  capable  of  accepting  values  of  y  from  30° 
to  145°  ;  then  a  may  vary  from  -55°  to  +60°  . 

Data  reduction  of  signals  for  a  typical  155  mm  M549  projectile 
in  free  flight,  using  a  yaw  sonde,  are  then  plotted  against  the  time 
required  to  study  the  dynamic  behavior  throughout  its  trajectory. 

Figures  4  and  5  show,  respectively,  the  spin  (PHI  DOT)  and  the  angular 
motion  (Sigma  N)  vs  time  (ref  10) . 

Time  Computation 

The  Greenwich  Mean  Time  (GMT)  of  a  firing  is  found  by  noting  the 
Local  Mean  Time  (LMT)  at  the  instant  of  firing  and  converting  as  follows: 

GMT  =  LMT  +  IZREF  (West  Longitude)  (2a) 

GMT  =  LMT  -  IZREF  (East  Longitude)  (2b) 

Values  of  IZREF  for  each  state  in  the  U.S.  are  given  in  table  1. 

For  locations  outside  the  U.S.,  the  value  of  IZREF  can  be  found  in 
reference  8.  As  an  example,  the  GMT  at  Yuma  Proving  Ground,  Arizona, 
for  LMT  equals  0830  is  found  to  be  1530  since  IZREF  equals  0700  according 
to  table  1. 

Celestial  Angle  Relations 

Figure  1  shows  the  angular  position  of  the  sun  and  the  observer 
with  reference  to  the  earth,  along  with  the  symbols  and  sign  conventions. 
The  positional  system  employed  here  (ref  12)  is  one  which  easily  defines 
the  position  of  an  observer  cn  the  earth  by  latitude  and  longitude. 

A  point  on  the  earth  directly  beneath  the  sun,  its  subpoint  (SS,  fig  1), 
could  just  as  easily  be  defined  by  its  latitude  and  longitude.  The 
usual  point  of  reference  for  longitude  (fig  1)  is  the  Greenwich  Meridian 
(GM)  0°  longitude.  The  longitude  component  of  the  position  of  the  sun's 
subpoint,  SS,  is  usually  expressed  as  a  hour  angle  either  from  the  Green¬ 
wich  Meridian  (the  Greenwich  Hour  Angle,  GHA)  or  from  the  observer's 
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meridan  (the  Local  Hour  Angle,  LHA) .  Longitude  is  measured  east  or 
west  from  0°  through  180° ,  but  hour  angles  are  always  measured  west¬ 
ward,  either  from  GM  or  the  observer* s  meridian  from  0°  through  360°. 

The  Local  Hour  Angle  is  computed  as  follows: 

LHA  =  GHA  -  West  Longitude  (3a) 

LHA  =  GHA  +  East  Longitude  (3b) 

Many  angular  measurements,  such  as  those  shown  in  figure  1,  are 
referenced  to  the  equator  (e.g.,  the  latitude  of  the  observer  and  the 
declination  (Dec)  of  the  sun)  or  to  the  celestial  horizon  for  the  ele¬ 
vation  of  the  sun.  The  declination  of  the  sun  is  also  the  latitude 
of  its  subpoint.  Latitude  is  measured  north  or  south  of  the  equator, 

0°  through  90°  at  the  poles.  Declination  is  limited  to  approximately 
23^  travel  north  or  south  of  the  equator  (due  to  the  orbit  of  the 
earth  about  the  sun  and  the  tilt  of  the  earth’s  north-south  axis). 

Finally,  we  need  to  define  the  angular  direction  in  which  to 
look  for  the  sun.  Figure  3  shows  that  the  azimuth  of  the  sun,ZN, 
from  an  observer’s  position,  is  measured  clockwise  from  north  to  a 
line  joining  the  observer’s  position  and  the  subpoint  of  the  sun. 

Determination  of  GHA  and  Declination  of  the  Sun 

Two  methods  of  finding  the  GHA  and  Dec  of  the  sun  are  presented. 

The  first  method  requires  a  simple  reference  to  the  Air  Nautical 
Almanac  for  a  specific  date  and  GMT.  [Table  2  is  extracted  from  the 
almanac  (ref  8).]  As  an  example,  for  28  Oct  1978  at  1530  GMT,  the 
almanac  states  that  the  GHA  equals  56°  32.5'  (56.542°)  and  the  Dec 
=  S  13°  08.4'  (-13.14°).  The  second  method  (ref  11)  employs  the  use  of 
table  3  and  is  explained  in  the  table.  It  is  an  emergency  method  of 
computing  the  sun’s  GHA  and  Dec  if  the  almanac  is  not  available.  At 
first  glance,  it  seems  cumbersome;  but,  with  practice,  the  computation 
is  accomplished  in  a  minute's  time.  Using  the  same  example  as  above  and 
following  the  explanation  in  table  3,  we  would  compute  GHA  equals  56°  33’ 
(56.55°)  and  Dec  equals  S  13°  09’  (-13.15°).  This  method  has  been 
programmed  into  the  computer  program,  SAZEL.  User  input  to  the  pro¬ 
gram  is  explained  in  appendix  A  and  a  listing  of  SAZEL  is  provided 
in  appendix  B. 

Determination  of  the  Sun's  LHA,  Elevation,  and  Azimuth 

Once  the  GHA  and  Dec  of  the  sun  have  been  found,  based  on  the 
date  and  GMT,  the  azimuth  (ZN)  and  elevation  (HC)  of  the  sun  can  be 
found  if  we  know  the  latitude  and  longitude  of  the  observer.  Through 
the  use  of  the  Sight  Reduction  Tables  for  Air  Navigation  (ref  11), 
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we  can  find  the  azimuth  and  elevation  of  the  sun.  (Table  4  is  an  ex¬ 
tract  from  reference  11  to  use  as  an  example.)  Using  the  GHA  and  Dec 
found  earlier,  and  placing  the  observer  at  Yuma  Proving  Ground,  lati¬ 
tude  equals  32.8811°  and  longitude  equals  114.3028°,  we  proceed  to 
find  the  LHA  using  equation  3a,  thus 

LHA  =  56.5500°  -  114.3028° 

LHA  =  302,2472?  (360°  has  been  added  to  make  LHA  positive) 

Entering  table  4  with  the  LHA,  Dec,  and  latitude,  we  find  the  azimuth 
and  elevation  of  the  sun  by  interpolation;  thus,  HO  =  18° 11'  and 
Z  =  ZN  =  120° .  The  computer  program,  SAZEL,  computes  HC  and  ZN  using 
rudimentary  principles  of  trigonometry. 

Six  Degrees-of-Freedom  Simulation  of  Free  Flight 
Missile  Using  SAZEL 

For  a  simulated  trajectory  of  a  free  flight  missile,  using  a 
modified  six  degrees-of-freedom  (6-DOF)  trajectory  program  (ref  13)  and 
SAZEL,  the  angular  motion,  a,  can  be  plotted  against  time  as  shown  in 
figure  6.  Figure  6  represents  a  simulation  of  the  155  mm  M549  projec¬ 
tile  which  was  fired  at  Yuma  Proving  Ground  (ref  10)  compared  to  the 
actual  recorded  angular  motion  shown  in  figure  4.  The  azimuth  and 
elevation  of  the  sun  can  be  represented  by  a  vector  S,  as  shown  in 
figure  3.  The  vector  S  can  further  be  represented  by  unit  vector  com¬ 
ponents  SI,  S2,  and  S3.  The  attitude  of  the  missile,  as  contained  in 
reference  13,  is  represented  by  the  vector  P.  Vector  P  is  also  repre¬ 
sented  by  its  unit  vector  components  PI,  P2,  and  P3.  The  dot  product 
of  these  unit  vectors  enables  us  to  find  y  and  finally,  a,  which  we 
are  ultimately  interested  in  plotting  against  time.  SI,  S2,  and  S3 
in  the  earth  frame  are  as  follows: 


SI  =  cos  HC  cos  ZN  (4a) 


S2  =  sin 

HC 

(4b) 

S3  =  cos 

HC  sin  ZN 

(4c) 

The  angle  y  is  then  computed  as  follows: 

y  =  arc  cos  (51P1  +  S2P2  +  S3P3) 
and  ultimately  we  find  a  from  equation  1. 
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Accuracy 


The  calibration  of  the  yaw  sonde  is  such  that  system  acquisition 
has  been  reported  to  be  ±  %°  .  The  method  of  computation  using  SAZEL 
is  reported  in  the  Air  Nautical  Almanac  to  be  ±  1'  (0.017°)  for  GHA 
and  ±  2'  (0.033°)  for  Dec  of  the  sun.  Small  errors  in  the  computation 
of  azimuth  and  elevation  would  be  inherent  due  to  GHA  and  Dec  from  the 
above  computation  and  changing  missile  position.  The  effect  of  changing 
latitude  and  longitude  after  launch  on  the  computation  of  GHA,  LHA,  HC, 
and  ZN,  can  be  minimized  by  using  the  latitude  and  longitude  at  the 
midpoint  of  a  trajectory.  Another  error  is  introduced  through  the 
necessary  assumption  of  unchanging  altitude  from  ground  level.  Finally, 
an  error  which  affects  both  computation  of  sun  elevation  and  data  acqui¬ 
sition  is  atmospheric  refraction.  Refraction  is  a  complicated  function 
of  height  above  the  earth  and  apparent  elevation  of  the  sun.  The  re¬ 
fraction  is  further  dependent  on  the  temperature  of  the  atmosphere  at 
the  height  of  observation.  At  any  given  apparent  sun  elevation,  the 
refraction  decreases  with  increasing  altitude.  Refraction  is  further 
decreased  if  the  atmosphere  is  warmer  than  normal  at  the  height  of  ob¬ 
servation.  For  the  flight  of  most  missiles,  the  refraction  encountered 
should  be  no  greater  than  -20'  (0.333°)  under  the  worst  conditions-- 
sea  level,  colder  than  normal  atmosphere,  and  sun  on  the  horizon. 

The  overall  error  for  ranges  under  30km  and  altitudes  from  ground 
level  to  30  km,  using  the  midpoint  trajectory  latitude  and  longitude, 
is  estimated  to  be  ±  ^  for  any  simulation. 

RESULTS 


The  reader  can  compare  the  manual  method  of  computation  for  LHA, 

GHA,  AZ  and  HC  to  that  derived  from  SAZEL.  Very  close  agreement  is 
shown  between  the  results  obtained  from  the  SAZEL  program  (app  A)  and 
those  extracted  from  the  Air  Almanac  and  Sightreading  Tables,  as  follows: 

Source  LHA  GHA  DEC  AZ  HC 

(computed) 

Air  Almanac  56.5417°  -13,1400° 


Sightreading  Tables  302.2472° 


56.5500°  -13.1500° 


Sight  Reduction  302.2472° 

Method 


120.0000°  18.1833° 


SAZEL 


56.5333° 


-13.1392°  119.8569°  18.2251° 


5 


Comparing  the  actual  angular  motion  in  figure  4  to  the  simulated 
angular  motion  in  figure  6,  close  agreement  is  obvious.  Note  that  the 
actual  flight  would  result  in  Sigma  N  from  the  apparent  sun  due  to  re¬ 
fraction.  A  suitable  aerodynamic  coefficient  package  and  the  test  data 
(record  of  firing  data  for  the  actual  flight)  are  necessary  to  obtain 
accurate  results. 


CONCLUSIONS 

SAZEL,  a  computer  program,  has  been  developed  which  determines 
the  azimuth  and  elevation  angles  of  the  sun  quickly  and  accurately, 
without  using  tabulated  sun  data  available  in  the  published  literature. 

In  comparison  with  the  yaw-sonde  method  (which  does  rely  on  previously 
published  data) ,  SAZEL  has  been  demonstrated  to  be  a  viable  tool  in 
simulating  the  angular  motion  of  missiles  during  free  flight. 

RECOMMENDATIONS 

1.  Update,  with  SAZEL,  the  ARRADCOM  computer  program  for  determ¬ 
ining  appropriate  times  (windows)  for  firing  projectiles  equipped  with 
yaw  sondes.  (The  trial  version  is  now  working.) 

2.  Complete  the  6-DOF  trajectory  program  using  SAZEL  which  allows 
computation  of  projectile  angular  motion  motion  (yaw  sonde  sun  angle.  Sigma  N) 
vs  time  for  plotting  on  the  ARRADCOM  graphics  terminal.  (The  trial  version 

is  now  working.) 
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Table  1.  IZREF  conversion  from  local  time  to  Greenwich  time 


STATE 

IZREF 

Alabama^- 

0600 

Alaska}  east  of  long. 

W.  137° 

0800 

long.  W.  137°  to  W. 

141° 

0900 

long.  W.  141°  to  W. 

161° 

1000 

long.  W.  161°  to  W. 

172° 

1100 

Aleutian  Islands 

1100 

Arizona 

0700 

Arkansas1 

0600 

California1 

0800 

Colorado1 

0700 

Connecticut1 

0500 

Delaware1 

0500 

District  of  Columbia1 
Florida1*2 

0500 

0500 

Georgia1 

0500 

Hawaii 

1000 

Idahc}*2 

0700 

Illinois1 

0600 

Indiana ^ 

0500 

Kansas^  * 2 

0600 

0600 

Kentucky*- » 2 

0500 

Louisiana1 

0600 

Maine1 

0500 

Maryland1 

0500 

Massachusetts1 

0500 

Michigan 

0500 

Minnesota1 

0600 

Mississippi1 

Missouri1 

0600 

0600 

Montana1 

0700 

Nebraska^  *  2 

0600 

Nevada1 

0800 

New  Hampshire1 

0500 

New  Jersey1 

New  Mexico1 

New  York1 

North  Carolina1 

North  Dakota*"*^ 

Ohio1 

Oklahom^1 

Oregorf  * 

0500 

0700 

0500 

0500 

0600 

0500 

0600 

0800 

Pennsylvania1 

0500 

Rhode  Island1 

0500 

South  Carolina1 

0500 

South  Dakota}  eastern 

part 

0600 

western 

part 

0700 

Tennesse^-*2 

Texas1 

0600 

0600 

0700 

Utah*-* 2 

Vermont^ 

Virginia1 

Washington,  D.C.1 
Washington1 

0500 

0500 

0500 

0300 

West  Virginia1 
Wisconsin1 

0500 

G600 

Wyoming1 

0700 

Summer  (daylight-saving)  time,  one  hour  fast  on  tl\e  time  given  (-100)  is 
kept  in  these  states  from  the  last  Sunday  in  April  to  the  last  Sunday  in 
October,  changing  at  02n  OO111  local  clock  time. 

2 

This  applies  to  the  greater  portion  of  the  state. 
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Table  2.  Extract  from  Air  Nautical  Almanac 


602  ftwmo})  GREENWICH  P.  M.  1978  OCTOBER  28  (SATURDAY) 


3MT 

'  o 

GHA 

Oec. 

ARIES 

GHA  °° 

VENUS-4.0 

GHA  Dec. 

JUPITER- 1.7 

GHA  Oec. 

SATURN 

GHA 

I  1.1 

Dec.  ■ 

9  MOON 

GHA  Oec. 

ZOO 

4  02.3 

S13 

05.5 

216  2E.6 

350  05  S24  04 

86  11  N18  41 

53  18  N 

8  52 

41  48  N 

2  05 

10 

6  32.3 

05.7 

218  59. C 

352  35 

88  41 

55  48 

44  13 

03 

20 

9  02.3 

05.8 

221  29.5 

355  06 

91  12 

58  19 

46  38 

02 

30 

11  32.3 

05.9 

223  59.9 

357  36  • 

93  42  • 

60  49  • 

49  04 

2  00 

40 

14  02.4 

06.1 

226  3C.3 

0  07 

96  12 

63  19 

51  29 

1  59 

50 

16  32.4 

06.2 

229  0C  7 

2  37 

98  43 

65  53 

53  55 

57 

00 

19  02.4 

S13 

06.4 

231  3Ll 

5  08  S24  04 

101  13  N18  41 

68  20  N 

8  52 

56  20  N 

1  55 

10 

21  32.4 

06.5 

234  0L5 

7  38 

103  43 

70  51 

58  45 

54 

20 

24  02.4 

06.6 

236  31.9 

10  09 

106  14 

73  21 

61  11 

52 

30 

26  32.4 

068 

239  02.3 

12  40  • 

108  44  • 

75  51  • 

63  36 

•  50 

40 

29  02.4 

06.9 

241  32.7 

15  10 

111  15 

78  22 

66  02 

49 

50 

31  32.4 

07.1 

244  03.2 

17  41 

113  45 

80  52 

68  27 

47 

100 

34  02.4 

S13 

07.2 

246  33  6 

20  11  S24  03 

116  15  N18  41 

83  23  N 

8  52 

70  52  N 

1  45 

10 

36  32.4 

07  3 

249  04  0 

22  42 

118  46 

85  53 

73  18 

44 

20 

39  02.4 

07.5 

251  344 

25  12 

121  16 

88  23 

75  43 

42 

30 

41  32.4 

07.6 

254  04  8 

27  43 

123  46  • 

90  54  • 

78  09 

•  40 

40 

44  02  5 

07  7 

256  35  2 

30  14 

126  17 

93  24 

80  34 

39 

50 

46  32.5 

07  9 

259  05  6 

32  44 

128  47 

95  54 

83  00 

37 

00 

49  02.5 

S13 

08.0 

261  36  0 

35  15  S24  03 

131  18  N18  41 

98  25  N 

8  52 

85  25  N 

1  36 

10 

51  32.5 

08.2 

264  064 

37  45 

133  48 

100  55 

87  50 

34 

20 

54  02.5 

08  1 

266  36  8 

40  16 

136  18 

103  25 

90  16 

32 

20 

56  32.5 

Q£},4 

?69  07  3 

42  47  • 

138  49  • 

105  56  • 

92  41 

31 

40 

59  02,5 

0&.6 

271  37  7 

45  17 

141  19 

108  26 

95  07 

29 

50 

61  32.5 

087 

274  08  1 

47  48 

143  49 

110  57 

97  32 

27 

00 

64  02.5 

S13 

08.9 

276  38.5 

50  18  S24  02 

146  20  N18  41 

113  27  N 

8  52 

99  57  N 

1  26 

10 

66  32  5 

09.0 

279  08.9 

52  49 

148  50 

115  57 

102  23 

24 

20 

69  02.5 

09.1 

281  39.3 

55  19 

151  21 

118  28 

104  48 

22 

30 

71  32.5 

09.3 

284  09.7 

57  50  • 

153  51  • 

120  58  • 

107  14 

•  21 

40 

74  02.6 

09.4 

286  40 J 

60  21 

156  21 

123  28 

109  39 

19 

50 

76  32.6 

09.6 

289  103 

62  51 

158  52 

125  59 

112  04 

17 

00 

79  02.6 

S13 

09.7 

291  41. D 

65  22  524  02 

161  22  N18  41 

128  29  N 

8  52 

114  30  N 

1  16 

10 

81  32.6 

09,8 

294  114 

67  52 

163  52 

131  00 

116  55 

14 

20 

84  02.6 

10.0 

296  41.3 

70  23 

166  23 

133  30 

119  20 

12 

30 

86  32.6 

10.1 

299  12.2 

72  53  • 

168  53  • 

136  00  • 

121  46 

11 

40 

89  02.6 

10  3 

301  423 

75  24 

171  24 

138  31 

124  11 

09 

50 

91  32,6 

10.4 

304  13.3 

77  55 

173  54 

141  01 

126  37 

08 

00 

94  02.6 

S13 

10.5 

306  43.4 

80  25  S24  01 

176  24  N18  41 

143  32  ^ 

8  52 

129  02  N 

1  06 

10 

96  32.6 

10.7 

309  13. i 

82  56 

178  55 

146  02 

131  27 

04 

20 

99  02.6 

10  8 

311  44.2 

85  26 

181  25 

148  32 

133  53 

03 

30 

101  32.6 

11.0 

314  14.’ 

87  57  • 

183  55  •  • 

151  03  • 

136  18 

1  01 

40 

104  02.6 

11  1 

316  45.. 

90  27 

186  26 

153  33 

138  44 

0  59 

50 

106  32.7 

11.2 

319  15.  > 

92  58 

188  56 

156  03 

141  09 

58 

00 

109  02.7 

S13 

11.4 

321  45.* 

95  29  S24  00 

191  27  N18  41 

158  34  N 

8  52 

143  34  N 

0  56 

10 

111  32.7 

11.5 

324  163 

97  59 

193  57 

161  04 

146  00 

54 

20 

114  02.7 

11  7 

326  46.“ 

100  30 

196  27 

163  34 

148  25 

53 

30 

116  32.7 

11.8 

329  17.1 

103  00  • 

198  58  • 

166  05  ■ 

150  51  • 

51 

40 

119  02.7 

11  9 

331  47.1 

105  31 

201  28 

168  35 

153  16 

49 

50 

121  32,7 

12.1 

334  17,* 

108  02 

203  58 

171  06 

155  41 

48 

00 

124  02  7 

S13 

12.2 

336  48.- 

110  32  524  00 

206  29  N18  41 

173  36  N 

8  52 

158  07  N 

0  46 

10 

126  32.7 

123 

339  18.6 

113  03 

208  59 

176  06 

160  32 

44 

20 

129  02  7 

12.5 

341  49.1 

115  33 

211  30 

178  37 

162  57 

43 

30 

131  32  7 

12.6 

344.  19. e 

118  04  • 

214  00  • 

181  07 

165  23  • 

41 

40 

134  02  7 

12.8 

346  50.C 

120  34 

216  30 

183  37 

167  48 

39 

50 

136  32  7 

12.9 

349  20a 

123  05 

219  01 

186  08 

170  14 

38 

00 

139  02.8 

S13 

13.0 

351  50. E 

125  36  S23  59 

221  31  N18  41 

188  38  N 

3  52 

172  39  N 

0  36 

10 

141  32  8 

13.2 

354  21.2 

128  06 

224  01 

191  09 

175  04 

34 

20 

144  02.8 

13.3 

356  51.6 

130  37 

226  32 

193  39 

177  30 

33 

30 

146  32.8 

13.5 

355  22.C 

133  07  • 

229  02  • 

196  09  • 

179  55  • 

31 

40 

149  02.8 

13.6 

1  52.5 

135  38 

231  33 

198  4C 

182  20 

29 

50 

151  32.8 

13.7 

4  22.9 

138  09 

234  03 

201  10 

184  46 

28 

00 

154  02.8 

S13 

13  9 

6  53  3 

140  39  S23  59 

236  33  N18  41 

203  41  N 

8  52 

187  11  N 

0  26 

10 

156  32.8 

14.0 

9  23.7 

143  10 

239  04 

206  11 

189  37 

24 

20 

159  02.8 

14.2 

11  54.1 

145  40 

241  34 

208  41 

192  02 

23 

30 

161  32.8 

14.3 

14  24.5 

148  11  • 

244  04  • 

211  12  • 

194  27  • 

21 

40 

164  02.8 

14.4 

16  54  9 

150  41 

246  35 

213  42 

196  53 

19 

50 

166  32.8 

14.6 

19  25.3 

153  12 

249  05 

216  12 

199  18 

18 

00 

169  02,8 

S13 

14.7 

21  55.7 

155  43  S23  58 

251  36  N18  41 

218  43  M 

8  51 

201  43  N 

0  16 

10 

171  32,9 

14.9 

24  26.1 

158  13 

254  06 

221  13 

204  09 

14 

20 

174  02.9 

15.0 

26  56.6 

160  44 

256  36 

223  43 

206  34 

13 

30 

176  32  9 

15.1 

29  27.0 

163  14  • 

259  07  ■ 

226  14  • 

208  59  • 

11 

40 

179  02.9 

15.3 

31  57.4 

165  45 

261  37 

228  44 

211  25 

09 

50 

181  32  9 

15.4 

34  27.8 

168  15 

264  07 

231  15 

213  50 

08 

o$* 

15  00.0 

SO  00.8 

15  03,5  NO  00,5 

15  02.2  SO  00.1 

15  02.3  SO  C0.1 

14  32.3  SO  09.9 

lot. 

Moon- 

set 

fa». 

:  n 

» 

00 

h 

i  * 

;  72 

15  37, 

-02 

i  70 

15 

.  00 

(  68 

15  39 

+  03 

r  66 

15  32 

05 

64 

15  31 

06 

>  62 

15  30 

07 

!  60 

15  29 

09 

I  58 

15  28 

10 

;  56 

15  27 

11 

-  54 

15  27 

12 

52 

15  26 

12 

50 

15  26 

13 

45 

15  25 

15 

40 

15  24 

16 

35 

15  23 

17 

30 

15  22 

18 

20 

15  21 

20 

'  10 

15  20 

22 

0 

15  18 

23 

10 

15  17 

25 

20 

15  16 

26 

30 

15  15 

28 

35 

15  14 

29 

40 

15  13 

30 

45 

15  12 

31 

50 

15  10 

33 

52 

15  10 

■  34 

54 

15  09 

35 

56 

15  08 

35 

58 

15  08 

37 

60 

15  07 

38 

S 

Moon's  P  in  A. 


Sun  SD  16. '1 
Moon  SD  15' 
Age  264 
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Table  3.  GHA  and  declination  of  the  sun  for  the  years  1965-2000 
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12 


c 

o 

•r- 

4-3 

ca 

CD 


> 

<X5 


U 

•r- 

ro 

S- 

o 

M— 

10 

a; 

-Q 

ra 

4-3 

c 

o 

4- 3 

U 

“O 

Q) 

5- 

4- 3 

CD 

•r- 

tO 

E 

o 

5- 
M— 


(/) 

4-3 

a 

ro 

s- 

4-3 

X 

LU 


Q) 


-Q 

fO 


13 


Sun 


West 


North 

Pole 


South 


Sun 


Celestial 

Horizon 


Pole 
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Figure  2.  Definition  of  gamma  and  sigma. 
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Figure  4.  Actual  sigma  N  vs  time  of  a  155  mm  M549  projectile. 
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Figure  5,  Actual  phi  dot  vs  time  of  a  155  ram  M549  projectile. 
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Figure  6.  Simulated  sigma  N  vs  time  of  a  155  mm  M549  projectile. 


APPENDIX  A 


INPUT  AND  OUTPUT  FOR 
COMPUTER  PROGRAM  "SAZEL" 


Description  of  Input  and  Output  of  SAZEL 

To  find  the  azimuth  and  elevation  of  the  sun  using  the  computer 
program,  SAZEL,  as  listed  in  Appendix  B,  data  are  input  on  one  computer 
card  as  follows: 


INPUT 


Column  No. 

Item 

Format 

1-2 

NDAY,  numerical  day  of  month 

12 

10-11 

NMONTH,  numerical  month  of  year 

12 

20-23 

NYEAR,  numerical  year  A.D. 
(.1978-2000  only) 

14 

30-31 

NHOUR,  GMT  hours 

12 

32-33 

NMIN,  minutes 

12 

50-59 

WLONG,  longitude,  degrees 

F9.4 

60-69 

WLAT,  latitude,  degrees 

F9.4 

EXAMPLE 

The  example  used  in  table  3  and  4  is  chosen  for  comparative 
purposes. 
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October  28,  1978;  GMT=  1530 

Yuma  Proving  Ground;  Longitude=114,3028  deg,  Latitude=32 , 8811  deg 


Column  No. 

Card  Ttem 

1-2 

28 

10-11 

10 

20-23 

1978 

30-31 

15 

32-33 

30 

50-59 

114.3028 

60-69 

32.8811 

OUTPUT 


ODAY  OTIME 

CZ)  GHA  AZIMUTH  DECLINATION  ELEVATION 

302  1600 

56.5333  119.8569  -13.1392  18,2251 
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LISTING  OF  COMPUTER  PROGRAM  "SAZEL" 
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GLOSSARY 


Azimuth 

Celestial  horizon 
Declination 

Elevation 

Greenwich  Hour  Angle 

Greenwich  Meridian 

Greenwich  Mean  Time 

IZREF 

Latitude 

Longitude 

Local  Hour  Angle 


The  angle  measured  clockwise  from  true  north 
to  a  line  passing  through  the  observer  and 
subpoint  of  the  sun 

The  tangent  plane  to  the  earth  passing  through 
the  observer's  position 

The  angular  distance  to  the  sun  measured 
north  or  south  through  90°  from  the  equator 
to  the  subpoint  of  the  sun 

The  angular  distance  of  the  sun  above  the 
celestial  horizon 

The  angular  distance  measured  from  the 
Greenwich  meridian  westward  through  360° 
to  the  meridian  passing  through  the  subpoint 
of  the  sun 

The  prime  meridian  which  passes  through 
Greenwich,  England,  from  which  longitude  is 
measured  east  or  west 

Local  time  at  the  Greenwich  meridian  measured 
by  reference  to  the  mean  sun.  It  is  the  angle 
measured  along  the  equator  (and  converted  to 
time)  from  the  Greenwich  meridian  westward 
through  360°  to  the  meridian  passing  through 
the  subpoint  of  the  mean  sun 

A  time  conversion  to  express  local  mean  time 
as  Greenwich  mean  time 

Angular  distance  measured  north  or  south  of 
the  equator  along  a  meridian,  0°  through  90° 

The  angular  distance  east  or  west  of  the 
Greenwich  meridian,  measured  along  a  line  of  ’ 
parallel  from  0°  to  180° 


The  angular  distance  measured  from  the  ob¬ 
server's  meridian  westward  through  the  subpoint 
of  the  sun 
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Local  Mean  Time 


Meridian 


Parallel 


Subpoint 


Time 


Z  or  Zulu  time 


Local  time  at  the  observer's  meridian  measured 
by  reference  to  the  mean  sun.  It  is  the  angle 
measured  along  the  equator  (and  converted  to 
time)  from  the  observer's  meridian  westward 
through  360°  to  the  meridan  passing  through 
the  subpoint  of  the  mean  sun 

Imaginary  line  on  the  earth  connecting  points 
of  equal  longitude 

Imaginary  line  on  the  earth  connecting  points 
of  equal  latitude 


That  point  on  the  earth's  surface  directly 
beneath  a  celestial  body 

Usually  expressed  in  four  numerals  (0001  thru 
2400)  where  there  are  24  hours  in  a  day  and 
60  minutes  to  one  hour.  6:29  PM  would  be 
expressed  as  1829 

An  expression  indicating  Greenwich  mean  time 
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Dec 

DOF 

GHA 

GM 

GMT 

HC 

IZREF 

LHA 

LMT 

OM 

^1  >  »  ^3 

Si  ,  S2  ,  S3 

ss 

ZN 

Y 


LIST  OF  SYMBOLS 

Declination  of  the  sun,  deg 

Degrees -of- freedom 

Greenwich  Hour  Angle  of  the  sun,  deg 

Greenwich  Meridian,  0°  longitude 

Greenwich  Mean  Time,  hours 

Elevation  of  the  sun,  deg 

Time  conversion,  hours 

Local  Hour  Angle  of  the  sun,  deg 

Local  Mean  Time,  hours 

Observer's  meridian,  deg  of  longitude 

Directional  cosines  of  missile  position,  P,  as  used  in 
6-DOF  trajectory  program  (earth  frame) 

Directional  cosines  of  sun  position,  S,  (earth  frame) 

Sun's  subpoint  meridian,  deg 

True  azimuth  of  the  sun,  deg 

Yaw  sonde  acquisition  angle  between  the  sun  and  pro¬ 
jectile  axis,  deg 

Angle  between  the  sun  and  transverse  plane  perpendicular 
to  projectile  longitudinal  axis,  deg 
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SIGN  CONVENTION 


Dec 

North  + 

South  - 

HC 

Above  celestial  horizon  + 
Below  celestial  horizon  - 

UZREF 

East  longitude  position  - 
West  longitude  position  + 

Latitude 

North  + 

South  - 

Longitude 

East  + 

West  - 

Sun 

Forward  of  projectile  + 
Aft  of  projectile 

ZN 

Clockwise  + 
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